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Abstract. We discuss the link between the nonthermal 
X-radiation and TeV 7-ray emission from SN 1006, and 
study the capabilities of both electronic and nucleonic 
models for explanation of the TeV flux observed from the 
northeast rim of SN 1006. We show that the interpretation 
of the TeV radiation by the inverse Compton scattering 
of electrons on 2.7 K cosmic microwave background radia- 
tion is possible, however due to the escape of high energy 
electrons, the 7-ray emission should be significantly con- 
tributed not only from the rim, but also from the inner 
parts of the remnant. This implies an angular size of the 
TeV 7-ray emission larger than the size of the nonthermal 
X-radiation. In this scenario the magnetic field in the rim 
should not exceed 10 /iG. Then, in order to allow accel- 
eration of particles well beyond 10 TeV the shock speed 
should be high, > 3000 km/s. The latter condition gives 
preference to a large distance to the source, c? ~ 2 kpc or 
so which is in the limits of distances currently discussed in 
the literature. On the other hand, a larger magnetic field 
of order 100/iG and smaller shock speeds (and therefore 
a small distance to the source of about 1 kpc) are not 
excluded. In that case the observed TeV radiation can be 
explained by shock accelerated protons in the rim through 
production and subsequent decay of Tr'^-mesons. Contrary 
to IC radiation, the 7r-decay 7-ray source coincides essen- 
tially with the rim, and therefore it could be recognized 
by a relatively compact angular width. Both the electronic 
and nucleonic models require high efficiency of particle ac- 
celeration close to the Bohm limit, and large total energy 
in accelerated particles at the level of respectively 1% and 
10% of the kinetic energy of explosion of SN 1006. We 
discuss observational possibilities to distinguish between 
electronic and nucleonic origins of 7-radiation. 
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1. Introduction 

The flux of cosmic rays is described by a smooth single 
power-law spectrum which extends up to the so-called 
'knee' around 10^^ eV. This can be interpreted as an ev- 
idence that the bulk of galactic cosmic rays is produced 
by a single source population. Since the early 60's the 
supernova remnants (SNRs) are believed to be the most 
probable sites of acceleration of galactic cosmic rays (see 
e.g. Ginzburg & Syrovatskii 1964), the basic argument be- 
ing that the galactic Supernovae (SNe), and the resulting 
supernova remnants (SNRs) are almost the only known 
potential sources with available (kinetic) energy needed to 
provide the observed luminosity of the Galaxy in cosmic 
rays, LcR > 10*°erg/s. Moreover, it has been shown (see 



e.g. Dru ry |1983| , Blandford & Eichl er [1987| , Berezhko k 
Krymski 198^ Jones & Ellison 1991) that a viable mech- 
anism - diffusive shock acceleration - can effectively oper- 
ate in relatively young SNRs accelerating particles to 100 



TeV (Lagage & Cesarsky 1983) or perhaps to higher ener- 
gies (Volk & Biermann 1988). The strong shocks in SNRs 
provide not only very effective, 10 per cent or more, con- 
version of the total SN explosion energy into the acceler- 
ated particles but also can naturally explain hard (cx E^^ 
with r ~ 2.0—2.1) source spectrum which follows from the 
standard 'leaky box' propagation model of galactic cosmic 
rays (see e.g. Gaisser 199C| ). 

The nonthermal (synchrotron) radio emission observed 
from shell-type SNRs is an unambiguous indicator of ac- 
celeration of GeV electrons there. Possible association of 
some of the 7-ray 'hot spots' detected by COS B (Bhat 
et al. |1985| ; Pollock |l985t Wo lfenda le & Zhang |l994i) and 
EGRET (Sturner & Dermer |1994| ; Esposito et al. ^996D 
with galactic SNRs could be another evidence for accel- 
eration of electrons and possibly also protons in SNRs 
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(Blandford & Cowie |1982|, S turner et al. |1997t Gaisser et 
al. 1998| ; Baring et al. 1999). However these observations 
concern only the low energy, typically i? < 10 GeV parti- 
cles, therefore they cannot yet conclusively argue in favor 
of SNRs as suppliers of the whole spectrum of galactic 
cosmic rays. An unbiased proof of this hypothesis can be 
provided only by detection of TeV 7-radiation produced 
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at interactions of accelerated protons and nuclei with am- 
bient gas through production and subsequent decay of 
7r°-niesons (Drury et al. |1994|; Naito & Takahara [1994D 



and/or by accelerated electrons upscattering the 2.7 K 
MBR (Mastichiadis |1996D . Therefore it is difficult to over- 
estimate the significance of the discovery of TeV radiation 
from the historical supernova SN 1006 reported by the 



CANGAROO collaboration (Tanimori et al. 1998) 



Actually the first, although rather circumstantial evi- 
dence for the presence of ultrarelativistic particles in SN 
1006 was found two decades ago by means of X-ray obser- 
vations. It was suggested that the power-law X-ray spec- 
trum of SN 1006 observed by the Einstein (Becker et al. 
1980| ) and EXOSAT (Jones & Pye |1989D satellites could 



be explained by synchrotron radiation of electrons ac 
celerated in the SNR shell (Reynolds & Chevalier 1981, 
Amosov et al. 1994). The recent detection of a spatially 



resolved power-law component of hard X-radiation from 
the edges of the remnant of SN 1006 by ASCA (Koyama 
et al. |1995D and ROSAT (Willingale et al. |1996D strongly 
supports the synchrotron origin of the radiation that im- 
plies acceleration of electrons to energies ~ 100 TeV. It 
should be noted that a featureless, power-law spectrum of 
X-radiation could be produced by a plas ma w ith multi- 
temperature structure (Hamilton et al. 1983 ), or by a 
heavy-element-dominated plasma with possible deviation 
of the electron distribution from the standard Maxwellian 



spectrum (Laming 1998 ) . Therefore the very fact of detec- 
tion of a power-law component of X-radiation does not yet 
automatically imply a nonthermal origin of the radiation. 
However in the case of SN 1006 the thermal interpretation 
of the continuous X-radiation faces serious difficulties, and 
the preference therefore is given to the nonthermal syn- 
chrotron hypothesis (Koyama et al. 1995; Laming 1998). 

Motivated by this fact, several theorists (Pohl |1996| ; 
Mastichiadis & de Jager |1996| ; Yoshida & Yanagita [l997| 
predicted detectable fluxes of inverse Compton (IC) TeV 
emission. Remarkably, very soon this object has been de- 
tected as a TeV emitter (Tanimori et al. 1998| ). Both the 
reported flux and the spatial position of the TeV source 
centered on the NE rim are close to the predictions. This 
circumstance strengthened the belief in the electronic ori- 
gin of TeV radiation. The 7r°-decay contribution to the 
observed flux is widely considered as less important. 

However, given the fact that SN 1006 is the only shell- 
type SNR reported as TeV 7-ray source, a careful exam- 
ination of different scenarios of 7-ray production is obvi- 
ously needed prior to exclude the nuclconic origin of the 



observed TeV radiation (Aharonian 1999). In particular, 
the contribution of the Tr^-decay 7-rays could be enhanced 
to the level of the observed fluxes if we assume a relatively 
small distance to the source of about ~ 1 kpc (Willingale 
et al. 1996), as well as take into account a significant com- 



pression of gas in the NE rim by the shock. It is impor- 
tant that this scenario does not limit the magnetic field 
in the acceleration region to values comparable with the 



interstellar B-field as required by the IC models of the ob- 
served TeV radiation. The magnetic field in the rim could 
be then as high a.s B ^ 100/iG which is favorable for accel- 
eration of particles up to ^ 10 TeV, taking into account 
the young age of SN 1006. Obviously in both scenarios the 
nonthermal X-ray emission observed from the NE rim is 
explained by synchrotron radiation of directly accelerated 
electrons. 

The current data do not allow us to give preference 
to the electronic or nucleonic origin of the observed TeV 
radiation. In this paper we study in detail the spectral 
and spatial characteristics of the IC and Tr^-decay 7-rays. 
In particular we show that while the size of 7r°-decay 7- 
ray emission is expected to be of the size of the rim, the 
size of the inverse Compton TeV source should be signif- 
icantly larger because of inevitable escape of very high 
energy electrons from the acceleration region, and their 
subsequent upscattering of the 2.7 K MBR outside of the 



2. Inverse Compton versus Tr^-decay 7-rays 

For a "standard" injection spectrum of shock-accelerated 
particles (protons and electrons) in a SNR, 



Q{E) cx exp (-E/Eq) , 



(1) 



with r « 2 and Eq « 100 TeV, the flux of 7r°-decay 7-rays 
at £: < 0.05£;o, 



F(> E) w 10"" ^(E/1 TeV)"^ ph/cm^s , 



(2) 



is determined by a simple scaling parameter, 

A = (lycR/lO^" erg) (n/1 cm^^) {d/l kpc)"^ ^ (3) 

Here Wcr is the total energy in relativistic protons, n 
is the ambient gas density, and d is the distance to the 
source. 

The peak luminosity of 7-rays is reached at the early 
Sedov phase (Berczhko & Volk |l997| ), i.e. typically 10^ - 
lO'' years after the SN explosion. At this stage the radius of 
the shell exceeds several parsecs. This results in a typical 
angular size of relatively close {d < 1 kpc) SNRs of about 
1°. The apparent conflict, from the point of view of source 
detect ability, between the angular size (0 oc 1/d) and the 
7-ray flux (_F!.y (X l/d^) significantly limits the number of 
SNRs which could be detected by the current most sen- 
sitive imaging atmospheric Cherenkov telescopes (lACT) , 
to the range of A > 0.1. Future large stereoscopic lACT 
arrays will extend this limit down to A ^ 0.01 (Aharonian 



Akerlof 1997) 



When extracting information about the accelerated 
protons one has to subtract a possible non-negligible 
"contamination" caused by parallelly accelerated electrons 
that upscatter the photons of 2.7 K MBR (which is the 
dominant target photon field for production of Te V pho- 
tons in most of SNRs; see e.g. Gaisser et al. |l998| ) up to 
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Fig. 1. The integral fluxes of 7r°-decay (heavy solid line) 
and IC (thin lines) 7-rays from a 10^ year old SNR. The 
tt" decay 7-rays are calculated for the scaling parameter 
A = \. The IC 7-ray fluxes are calculated for 3 differ- 
ent values of the magnetic field B = i (solid curve), 10 
(dashed curve), and 30 /iG (dot-dashed curve), assuming 
the electrons produce the same flux of synchrotron radia- 
tion Si, = 10 /^Jy. For both the protons and electrons we 
assume the same injection spectrum given in the form of 
Eq.(l) with F = 2 and -Bo = 100 TeV. 

7-ray energies. For production of IC 7-rays of TeV ener- 
gies multi-TeV electrons are needed. The same electrons 
produce also synchrotron UV/X-ray radiation. The typi- 
cal energies and E'x of the IC and synchrotron photons 
produced by an electron are related as (see e.g. Aharonian 
et al. ^997|) [| 

(4) 



~ 1.5(£;x/0.1kcV) (B/10AiG)"iTcV. 



The ratio of the synchrotron and IC fluxes /e = 
E'^F{E) — vSi, at these energies does not practically de- 
pend on the shape of the spectrum of parent relativistic 
electrons, but strongly depends on the magnetic field: 



fsy{Ex) 



0.1 {B/iOfxcy 



(5) 



For the X-ray spectrum with photon index ax ~ 2 the 
energy flux fx is almost energy-independent. Therefore fx 
at a typical energy of 1 keV could serve as a good indicator 
for the IC 7-ray fluxes expected at TeV energies, although 
for magnetic fields B < 100 /iG the energy of synchrotron 
photons relevant (produced by the same parent electrons) 
to ~ 1 TeV 7-rays is in the soft X-ray domain (see Eq. 4). 

The contribution of 7r°-decay 7-rays dominates over 
the contribution of the IC component when 



A > 0.1 (Sikev/lO^Jy) (B/IO AiG)" 



(6) 



where SikoV is the flux of nonthermal synchrotron ra- 
diation at 1 keV normalized to 10/iJy (the correspond- 
ing energy flux fx ~ 2.4 x 10"^^ erg/cm^s). Note that 
SikoV — 10 /iJy is a typical level of nonthermal X-ray 
fluxes reported recently from four shell-type SNRs - SN 



1006 (Koyama et al |1995| ), Cas A (Allen et al. [1997D , IC 
443 ( Keoha ne et al. |1997D and RXJ1713. 7-3946 (Koyama 
et al. \m% . 

In Fig. |l| we present integral fluxes of Tr'^-decay and in- 
verse Compton 7-rays from a SNR of age 10^ yr. The 7r°- 
decay 7-ray flux corresponds to the scaling factor A — 1. 
The IC fluxes are calculated by normalizing the syn- 
chrotron X-ray fluxes to SikeV = lOA'Jy for ambient B- 
fields 3/iG, 10/LtG, and 30/xG. For both electrons and pro- 
tons we assume continuous acceleration during 10'^ years 
with a time-independent source spectrum taken in the 
form of Eq.(l) with F = 2 and = 100 TeV. Note that 
for the normalizations used, the results presented in Fig. |^ 
only slightly depend on the source age, unless it is larger 
than the radiative cooling time of multi-TeV electrons, 
irad « 10^ (B/10AiG)-2(E^/ioTeV)-iyr. 

3. Synchrotron and IC components of nonthermal 
radiation of SN 1006 

The bulk (~ 75%) of the nonthermal X-ray component 
of SN 1006 is contributed by the northeast (NE) and the 
southwest (SW) rims (Koyama et al. 1995, Willingale et al. 
1996| ). The observations of SN 1006 by the CANGAROO 
collaboration (Tanimori et al. 1998) have revealed TeV 
radiation from the NE rim, but only an upper limit from 
the SW rim was reported. Therefore in this paper the 
nonthermal radiation of only NE rim is discussed. 

For the spectral characteristics of the nonthermal X- 
ray emission of NE rim we follow the procedure used by 



^ Note that this relation neglects the Klein-Nishina effect 
which becomes important at energies > 10 TeV. 



Mastichiadis & de Jager (1996). Namely we assume that 
the X-ray spectrum below 1.5 keV is described by a power- 
law with mean photon index — 2.2 as it follows from 
the combined Einstein/EXOSAT/ROSAT data. At higher 
energies up to 8 keV the spectrum reported by ASCA from 
NE-I-SW rims is significantly steeper, with a power-law 
photon index ax — 3. For the estimate of the absolute 
X-ray flux from NE rim we take into account that ~ 60% 
of the total emission of NE and SW rims is contributed 
by NE rim, as it follows from t he res pective count rates of 
ROSAT (see Willingale et al. |1996[ ). The resuhing range 
of the X-ray fluxes observed from the NE rim in the range 
from 0.1 keV to 8 keV is shown in Fig. |[ 

The radio emission of SN 1006 has a power-law spec- 
trum with an index a^ ~ 0.57 and intensity 30.8 Jy at 
408 MHz (Stephenson et al. |l977| ). In Fig. | one half of 
this flux is taken because the observed total radio flux 
is equally contributed from the northeast and southwest 
directions (see e.g. Reynolds & Gilmore 1986 ). 

In Fig. H we show also the range of TeV 7-ray fiuxes 
observed by CANGAROO from the direction of NE rim 
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of SN 1006 at energies above 1.7 TeV (Tanimori et al. 
1998| ), as well as the EGRET flux upper limit from SN 



1006 (Mori 1997) 



3.1. One zone model 

In this section we assume that synchrotron and IC 7- 
rays are produced by relativistic electrons confined in a 
spatially homogeneous rim. In this simplified case the IC 
fluxes shown in Fig 1 could serve as first approximation es- 
timates for the expected TeV 7-radiation. However, more 
reliable predictions of 7-ray fluxes require a careful deriva- 
tion of the spectrum of > 10 TeV electrons which is con- 
trolled by the spectrum of synchrotron X-rays. 
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Fig. 2. The synchrotron (heavy lines) and IC (thin lines) 
fluxes calculated for the homogeneous source without es- 
cape of accelerated particles (one-zone model) . Power-law 
injection spectrum of electrons with Fe = 2.15 is assumed 
in order to flt the radio data. The maximum energy Eq 
is determined from the condition given by Eq.(9) for 3 
different values of the magnetic held: B = 3 (dashed 
line), 5/LiG (solid line), and 8/xG (dot-dashed fine). 

In the synchrotron-inverse Compton models the spec- 
tral fit to the X-ray flux is crucial because the fluxes are 
produced by electrons in the region of the exponential 
cutoff Eq between 10 and 100 TeV. The spectral index 
F of accelerated electrons is derived from the radio data, 
F = 1 + 2ar ~ 2.15. Information about Eq in Eq. (1) is 
contained in the X-ray spectrum . In particular case of 
negligible energy losses, the spectrum of electrons N{Ec) 
repeats the injection spectrum, N{Ee) oc Q{Ec), there- 
fore in the (5-functional approximation for the synchrotron 
emissivity, the X-ray spectrum can be pre sente d in a sim- 
ple f orm (see e.g. Yoshida & Yanagita 1997 , Reynolds 
199^ ) 

FiE) oc ii;-(r+i)/2 exp[-(i?/i?^)i/2] , (7) 
where 

E,„ ~ 0.29hvr 



is the mean energy of synchrotron photons produced by an 
electron of energy Eq [vc is the characteristic synchrotron 
frequency ). 

From Eq. ^ follows that the spectrum of synchrotron 
radiation depends only on the product EqB^/"^. Numer- 
ical calculations using the accurate expression for syn- 
chrotron radiation generally confirm that approximation 
in the form of Eq.(7) holds at least until E ^ 10i?ni 
(Reynolds [T998| ). 

In Fig. 2. we show 3 spectra of synchrotron radiation 
calculated for the same product 



n = EB^/^ = 61.5 TeV /iG 



1/2 



(9) 



but for 3 different combinations of Eq and B: (1) B = 
SfiG, Eq = 35.5 TeV; (2) B = bfiG, Eq = 27.5 TeV; (3) 
B = S^iG, Eq = 21.3 TeV 2. In all calculations we use the 
same normalization to the radio fiuxes. Therefore different 
magnetic fields require different total energy in relativistic 
electrons. At the same time, because the target photon 
field (2.7 K MBR) for the IC scattering is fixed, there is a 
strong dependence of the IC 7-ray fluxes on the magnetic 
field. In particular, in the case of small magnetic fields 
B < 10/iG, when the radiative cooling time 



1.2 



lO'^ (E/IO TeV)"i(B/10 fiG)-^ yr 



(10) 



of 10 to 100 TeV electrons responsible for the observed 
X-rays and TeV 7-rays exceeds the age of the source, we 
have F^ oc ^-("'+1). Note that for B < 3/^G the radiative 
losses on 2.7 K MBR become dominant. 

Fig. ^ shows that the interpretation of the observed 
TeV radiation by IC mechanism in the framework of a 
simplified spatially homogeneous (one-zone) model is pos- 
sible only for an ambient magnetic field in a very nar- 
row range around 5/iG. Magnetic fields > 7/xG result in a 
strong reduction of the IC flux, while the assumption of 
low magnetic fields, B < 4/xG, leads to overproduction of 
7-rays. 

For such small magnetic fields a question arises 
whether the acceleration of electrons could be efficient 
enough in order to boost particles up to energies Eq — 
n/B^/^ ^ 25 — 30 TeV during the time of acceleration to 
limited by the age 10'^ yr of SN 1006. Indeed, in the model 
of diffusive shock acceleration the characteristic maximum 
energy Eq that a particle could achieve during the time Iq 



is determined (see e.g. Lagage and Cesarsky 1983) by the 
shock velocity Vs, the magnetic field B, and the so-called 
gyrofactor rj (the ratio of the particle mean free path to 
the gyroradius): 



Eq 



20 



B 



lOyuCy Vio^yr 



to 



0.02 (B/lO^iG) [Eq/IO TeV)2 keV (8) 



^ Note that our best fit value for the product II = E B^^^ = 
61.5TeV ^G^^'^ is b y a f actor of two less than obtained by 
Yoshida & Yanagita 1997| . The reason presumably is connected 
with different fluxes of X-rays used in calculations. 
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^2000 km/s 



(11) 



This estimate would be not significantly changed even if 
we take into account that in the past the shock speed 
was much higher than the present value Vg in Eq.(ll). 
Indeed, even assuming that the Sedov phase in SN 1006 
has started at ^ 1000 yr and applying the Sedov solu- 
tion for the shock speed in the past, v{t) = Vs{t/tQ)^^^^ 
for t > t^, the maximum energy of accelerated particles 
which is reached at f could be higher than it is given 
by Eq.(ll) only by a factor of (to/t*)^/^. Therefore, even 
for ~ O.lio ~ 100 yr, the maximum energy of the par- 
ticles in the past could be only by a factor of 1.5 higher 
than it follows from Eq.(ll). On the other hand Eq.(ll) 
does not include possible dissipative effects, such as the 
wave damping and the synchrotron losses of the electrons, 
therefore it can be considered as a good estimate for the 
maximum possible energy of accelerated particles. 

Eqs.(9) and (11) result in a lower limit on the strength 
of the magnetic field 



B> 12 



tn 



103 yr 



-2/3 



2000 km/s 



-4/3 



r;2/3;xG (12) 



The Sedov solution for the shock speed in SN 1006 
gives (Willingale et al. 1996| ) 



2000(d/lkpc)2km/s. 



(13) 



The shoc k spe ed 1400 ± 200km/c estimated by Will- 
ingale et al. ( |1996| ) from the ROSAT X-ray data implies a 
rather small distance to the source, d = 0.7±0.1 kpc. Win- 



kler & Long (1997) argued that the distance to the source 
should be larger, d = 1.8 ±0.3 kpc. At the same time, the 
shock spe ed Vs = 2600 ± 300km/c preferred by Winkler 
& Long ( 1997 ) is significantly lower than ~ 3600 km/c 



following from Eq.(12) at d = 1.8 kpc. Given the current 
controversy both in the estimated shock speed and the dis- 
tance, below we will discuss distances in the range from 
0.7 to 1.8 kpc, but using for all distances Eq.(13). 

Thus, if the shock speed does not significantly exceed 
2000 km/s which implies a distance ~ 1 kpc, even in the 
Bohm limit {rj = 1) the strength of the magnetic field 
should not be less than 10 /xG even if we assume that the 
effective particle acceleration has started at early stages 
of this SNR. This limitation could be softened by a factor 
of 2 for larger distances to the source (see below). 

Another way to soften the conflict between the ses- 
timate of the B-field given by Eq.(12) and the value of 
B < 7 fiG needed for explanation of TeV fluxes is pos- 
sible if we assume that the bulk of the observed 7-ray 
flux is produced outside of the shell. The highest energy 
electrons produced in the shell can effectively escape the 
acceleration region. This will give rise to an enhanced IC 
emission outside of the rim, namely in the interior regions 
of the remnant where the magnetic fleld could be as low 



as the typical interstellar B-field. Such a leakage of elec- 
trons from the rim is unavoidable because of diffusive and 
convective propagation of particles. This effect should be 
taken into account not only for correct derivation of the 
electron spectrum inside the rim, but also for calculations 
of the IC radiation outside of the rim. Below we show that 
in fact the contribution of this component (missed in pre- 
vious studies) is comparable or even can exceed TeV flux 
produced in the rim. 

3.2. Inverse Compton TeV emission outside of the rim 

The radio and X-ray images of SN 1006 indicate that the 
most probable site of particle acceleration is the thin shell 
of the remnant, in particular the NE rim. The width Ar of 
the NE rim in X-rays as measured by ROSAT does not ex- 
ceed 20% of the angular radius of the remnant of about 17 



arcminutes (Willingale et al. 1996) which corresponds to 
Tg = 4.9 dkpc pc at a distance to the source d — 1 d^pc kpc. 
In calculations we assume a constant Ar/rg = 0.2 ratio 
throughout the evolution of the remnant. 

Due to two principal mechanisms of propagation - dif- 
fusion and convection - the relativistic particles cannot be 
completely confined in the rim. The characteristic escape 
time is t^sc = {t^oI + ^dif 
is 



('^con + "^dif ) ^ ■ The convective escape time 



,(i) = Ar/u2 



(14) 



where U2 is the fiuid speed downstream of the shock (in 
it's rest frame), thus at present U2 = Vs/ p ^ 500 km/s for 
the compression ratio p ^ A. The width of the shell Ar ~ 
0.2rs oc t^/^ assuming a Sedov- type solution for the shock 
radius rg. Note that strictly speaking this solution is valid 
only at stages after the onset of the Sedov phase, i.e. for 
t* < i < to- However the energy distribution of electrons 
formed in the shell at present should be dominated by 
particles accelerated recently, t > 0.3to or so, unless one 
would assume that most of particles have been injected 
at very early stages. Therefore possible deviation of Ar{t) 
from the law oc t'^^^ at times t -C to does not significantly 
affect the spectrum of electrons at present. 
The diffusive escape time is 



Tdiff(t) 



(Ar)^ 
2D{E) 



(15) 



where D{E) is the diffusion coefficient in the shocked re- 
gion (the rim). In the theory of shock acceleration the 
diffusion coefficient is generally taken in the form 

D{E)^/-^, (16) 

where r^ = 3.3 • 10i5(£;/10 TeV)(B/10 ^iG)"! cm is the 
particle gyroradius, and 77 > 1 is the gyrofactor. The max- 
imum confinement (and therefore the maximum energy) 
of particles is achieved in the Bohm diffusion regime cor- 
responding to r] — 1. 
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The kinetic equation for energy distribution of rela- 
tivistic electrons N{E) in a source reads (see e.g. Ginzburg 
& Syrovatskii[l964|): 



dN 



d 



dt dE, 



[PN] 



N 



Q, 



(17) 



where Q = Q{Ec,t) is the electron injection rate, and 
P = —dEc/dt is the electron energy loss rate. The solution 
to this equation, including the time an d ene rgy dependent 
escape losses is (Atoyan & Aharonian 1999| ): 



NiE,,t) 



1 



P{E, 
exp 



P{Ct)Q{Ct,h) X 



da; 



ti '''esc (Ca; ? 2!') 



dti 



(18) 



The variable Q corresponds to the initial energy of an 
electron at instant ti which is cooled down to given energy 
Ec by the instant t, and is determined from the equation 



dE 



P{E) 



(19) 



In this paper we consider only overall (i.e. integrated 
over the volumes) fluxes from the NE rim (zone 1) and 
outside (zone 2). The electron distribution in the rim 
Ni{Ee, t) is found assuming injection of electrons at a con- 
stant rate during the age of SN 1006, with the spectrum 
Qi{Ee) in the form of Eq.(l). Although the injection rate 
of accelerated electrons is actually time-dependent, the 
main contribution to the total amount of particles in the 
rim comes from the acceleration at epochs from t ^ t^, 
to several times (t* corresponding to the onset of the 
Sedov phase) when the injection rate Qi{t) is almost time- 
independent (see Berezhko & Volk |1997| ). Since SN 1006 
most probably has already reached the Sedov phase (Volk 



1997) an assumption of time-independent Qi{E,t) is an 



adequate approximation (see below). At the same time we 
assume time-dependent escape, as described above (see 
Eq. 14 and 15), which have a stronger impact on the re- 
sults than the weak time dependence of the injection rate. 

Outside the rim the energy distribution N2{Ec,t) is 
found from the same Eq.(18) where the injection function 
is given by the rate of escape of electrons from the rim, 
Q2{Ec,t) = Nx/tssc- The region outside of the rim cor- 
responds to the interior of the remnant close to NE rim. 
Therefore in this region we ignore the losses due to escape. 

In Fig. H we show the synchrotron and IC fluxes from 
the NE rim and from the inner part of the remnant for the 
magnetic fleld in the rim Bi = 8 fxG and B2 = 2> /iG in the 
remnant. In order to show the impact of diffusive escape 
on the resulting radiation spectra, we considered 2 values 
for the parameter rj. The case of maximum confinement of 
particles in the rim, = 1, is shown by solid lines, while 
the dashed lines correspond to the case of 77 = 10. The 



77 — 1 (dashed) 
7) — 1 (solid) 



EGRET CANGAROO 




Log(E/eV) 

Fig. 3. The synchrotron and IC fluxes produced inside 
(heavy lines) and outside (thin lines) of the rim calculated 
within 2-zone model for two values of the gyrofactor: t] — 1 
(solid), and 10 (dashed). The magnetic fields inside and 
outside of the rim are taken as i?i = 8 fiG and ^2 = 3 /iG, 
respectively. The distance to the source is assumed 1 kpc. 



heavy and thin lines correspond to the emission produced 
inside and outside of the rim, respectively. 

For the parameters described above the observed X- 
ray spectrum of the rim is best explained for the expo- 
nential cutoff energy Eq = 30 TeV. The increase in the 
value of Eq compared with the best-fit value of 21.3 TeV 
used in Fig. 2 for the same magnetic field 8 /xG in the 
simple single-zone approach, is explained by diffusive es- 
cape of the electrons from the rim. For the magnetic field 
used, even in the Bohm limit the diffusive escape time of 
electrons at present is Tdif ~ 3500 yr/(E/10 TeV), there- 
fore for electrons with energy E ^ (30 — 100) TeV which 
are responsible for the X-ray fluxes, the diffusive escape 
time becomes comparable or even less than the age of the 
source. This results in a modiflcation (steepening) of the 
spectral shape at high energies. In order to compensate 
this effect, a higher Eg is to be assumed. Note that the 
convective escape time at present, T^on — 2000 yr, is also 
comparable with the age of the source. However the con- 
vective escape does not modify the energy spectrum of 
electrons. This is true, to some extent, also for radiative 
losses in the case of small magnetic fields B ^ 10 /iG, since 
the synchrotron cooling time then equals the source age 
only at E^ > 100 TeV. 

The electron escape from the rim leads to nonthermal 
X-ray production in a broader region (zone 2). Note that 
for the case of diffusion in Bohm limit (77 = 1) and the 
assumed magnetic field B2 = 3/iG, the X-ray fiux out- 
side of the NE rim makes ~ 10% of the flux produced 
in the rim. This implies that the magnetic fleld B2 can- 
not be larger than 5 /xG, otherwise the ratio of the fluxes 
produced outside and inside of the rim would exceed the 
value of about 1/3 as observed by ASCA (Koyama et al. 
19951) and ROSAT (Wilhngale et al. |l996| ). 
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While the synchrotron radiation outside of the rim can 
be suppressed assuming low magnetic field B2, the forma- 
tion of an extended IC radiation is unavoidable due to the 
uniform distribution of the principal target photon field 
for the Compton scattering - 2.7 K MBR. In particular 
even for 77 = 1 the contribution of 7-ray fluxes at > 1 TeV 
from the region outside of the rim becomes comparable 
with the flux from the rim. In the case of fast diffusion 
corresponding to rj = 10, the contribution from the ex- 
tended region outside of the rim, i.e. inside the remnant at 
r < O-Sts, well dominates the overall TeV flux of SN 1006 
(see Fig. I). 

The value of Eq = 30 TeV in Fig. 3 is by a factor 1.5 
larger than the estimate of the characteristic maximum 
energy of electrons following from Eq. (11) for the shock 
speed 2000 km/s corresponding to a distance to the source 
d = Ikpc (see Eq. 13). This conflict could be overcome 
assuming a higher speed of the shock, i.e. lager d. It is im- 
portant to note that the assumed distance to the source 
has a significant impact also on the timescale of the diffu- 
sive (but not the convective) escape of the electrons from 
the shell as Tdif cx Ar^ oc d^ (see Eq. 15). Therefore the 
relative contribution of the TeV emission produced in the 
interior of the remnant depends on the assumed distance 
to the source. This effect is seen in Fig. 4 where we show 
the results of calculations for two different distances cur- 
rently discussed in the literature, d = 0.7 kpc (WiUingale 
et al. |1996D and d = 1.8 kpc (Winkler & Long |l997D . For 
the distance d = 0.7 kpc we used Eq = 18 TeV which is 2 
times larger than it follows from the estimate of the max- 
imum energy by Eq. (11). However, even under such an 
extreme assumption the calculated synchrotron fluxes roll 
over too early, failing to reach the observed X-ray fluxes. 
Meanwhile for the distance d = l.Skpc the chosen value of 
Eq = 35 TeV is in agreement with Eq (11) even for ry = 2. 
We note that the contribution of the TeV flux produced 
outside of the rim is still very significant even for this large 
distance. 

The existence of the external IC component due to the 
escape of electrons from the rim results in a significant 
increase of the overall TeV radiation, and thus allows a 
larger magnetic field as compared with the value Bi = 
5 ± 1 fiG derived within the one-zone model (Section 3.1). 
In order to investigate the impact of this effect we have 
calculated the synchrotron and IC fluxes for two different 
combinations of Bi and 77 shown in In Fig. 5. In both cases 
we assume for the distance to the source d = l.Skpc. The 
maximum electron energy i?o = 32 TeV is calculated from 
Eq.(ll). We also assume low magnetic field outside of the 
rim, B2 = Bi/A, in order to avoid overproduction of X- 
radiation in this region. 

The combination -Brim = -Bi 5 /iG, t] — I provides 
a good fit to the X-ray data, but it leads to an over- 
production of the total TeV emission compared with the 
CANGAROO measurements by a factor of two. Given a 
possibility that the TeV emission from extended region 




Log(E/eV) 

Fig. 4. The synchrotron and IC fluxes produced inside 
(heavy lines) and outside (thin lines) of the rim calculated 
within 2-zone model for two distances to the source: d = 
1.8 kpc, and d — 0.7 kpc. The gyrofactor and maximum 
energy are: 77 = 2 and Eo — 35 TeV for d — 1.8 kpc (solid 
lines), and 77 = 1 and Eq — 18 TeV for d = 0.7 kpc (dashed 
lines). The magnetic fields inside and outside of the rim 
are taken as i?i = 6 /iG and B2 = 1.5 /iG, respectively. 



of the remnant (outside of the rim), which contributes 
half of the total flux, could be difficult to extract by the 
CANGAROO, this combination of parameters seems still 
acceptable. The assumption of Bi = 10 /iC and 77 = 2 
results in overproduction of X-rays by 50%. At the same 
time the predicted fluxes of TeV 7-rays from the rim are 
below the reported fluxes by a factor of 3. 

From Fig. 5 we may draw a conclusion that the mag- 
netic field in the rim should be within the limits 5 < i3i < 
10 /iG (preferably Bi ~ 6 — 8/xG), and 77 close to 1, e.g. 
the diffusion in the rim should take place essentially in 
the Bohm limit. In order to provide maximum electron en- 
ergy Eq of order of 30 TeV, the shock speed should exceed 
3000 km/s which implies a large distance to the source, 
d > 1.5 kpc. For this set of parameters, approximately half 
of the total TeV emission is contributed from the inner 
parts of the remnant, r < O.Sts. The narrow range of the 
required magnetic field in the rim allows rather accurate 
estimate of the total energy of relativistic electrons in SN 
1006. For example for the distance d = l.Skpc in Fig. 4 
the total electron energy in the rim is We — 2.9 x lO^^erg, 
whereas the energy in the electrons escaped the rim is 
Wo — 1.9 X lO'^^erg. It should be noted that the spectrum 
of electrons outside of the rim is enriched by multi-TeV 
particles due to the effect of energy dependent (oc E) dif- 
fusive escape. The total energy in electrons required to 
explain the X-ray and TeV 7-ray fluxes is about 5 x 10^* 
erg. Since the total energy of electrons depends on the B- 
field as Wo oc the total electron energy could be 

reduced only by a factor of (8/6)^'^ ~ 1.6 for the magnetic 
field B = 8 /iG. Therefore the inverse Compton origin of 
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Fig. 5. The synchrotron (top panel) and IC 7-ray (bottom 
panel) fluxes calculated in the framework of 2-zone model 
for two values of the magnetic field in the rim: Bi = 5 [iQ 
(solid lines) and B\ = 10 /xG (dashed lines). The mag- 
netic field outside of the rim is taken as Bi = Bi/A. The 
heavy lines show the fluxes from the NE rim, and thin 
lines correspond to the total fluxes. The distance to the 
source is assumed 1.8 kpc. The maximum electron energy 
Eq is calculated from Eq. (11) for the gyrofactor r/ = 1 for 
Bi = 5/iG, and ?7 = 2 for Bi = lO/xG. 

the detected TeV radiation requires about 1 percent of the 
total explosion energy of SN 1006 in relativistic electrons. 

The total energy of the magnetic field in the rem- 
nant with B2 = Bi/A ~ 2/iG for the distance 1.8 kpc 
is about 10^^ erg. A similar amount of magnetic field en- 
ergy is expected also in the rim since the amplification 
of the B-field there is compensated by a smaller volume 
of the rim. Thus, the inverse Compton origin of TeV 7- 
radiation implies that the conditions in SN 1006 are far 
from the equipartition between relativistic electrons and 
the magnetic field. 

In all calculations above we have assumed a station- 
ary injection rate of relativistic electrons which is a good 
approximation as far as the source is in an early Sedov 
phase. In order to demonstrate that, in Fig. 6 we show 
the spectra of synchrotron and IC fluxes calculated for 



two different assumptions: (i) a timc-indcpcndcnt injec- 
tion rate Q(t) = const, (ii) injection rate in the form 
Q{t) oc {l+t/tt,)~^ with a characteristic time (onset of the 
Sedov phase) = 100 yr. In both cases we have assumed 
acceleration spectrum given by Eq. (1) with Eq = 30 TeV. 
For the assumed magnetic field in the rim Bi = 7 and 
the current shock speed = 3600 km/s this value of Eq 
implies a gyrofactor ij = 1.5 (see Eq. 11). As it is seen 
from Fig. 6 the higher injection rate in the past does not 
affect significantly the fluxes of nonthermal radiation even 
for a very small values of = 100 yr. The main difference 
between the stationary and the time-dependent injection 
cases is reduced to somewhat higher IC 7-ray fluxes out- 
side of the rim. For more realistic values of > 300 yr 
the effect of time-dependent injection becomes practically 
negligible. 
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Fig. 6. The synchrotron (top panel) and IC 7-ray (bottom 
panel) fluxes calculated in the framework of 2-zone model 
for two assumptions about the injection rate of electrons: 
Qif) = const (solid lines), and Q{t) oc (1 -|- t/t»)~^ with 
U = 100 yr. In both cases B = 7/iG and Eq = 30 TeV. 
The distance to the source is 1.8 kpc. 

The above calculations show that the electronic origin 
of the observed TeV 7-ray emission meets difficulties for 
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small distances to the source, d ~ 1 kpc or less, since the 
latter implies a small shock speed, and therefore insuffi- 
cient acceleration rate in order to accelerate electrons to 
energies of about 30 TeV and beyond. In case of confir- 
mation of the recent claims about small distance to the 



source (Willingale et al. 1996), a possible solution to this 
difficulty could be an assumption of a more effective ac- 
celeration mechanism compared with the conventional dif- 
fusive shock acceleration, e.g. by invoking the mecha nism 
of perpendicular diffusion proposed by Jokipii 1987 . An- 
other principal possibility to increase the accelerate rate 
is connected with an assumption of a large magnetic field 
in the rim, Bi ^ lOfiG. Obviously in this scenario the 
TeV radiation principally cannot be explained by the IC 
mechanism. The only reasonable alternative could be the 
nucleonic origin of TeV radiation. 



4. TeV gamma rays of nucleonic origin 

Currently, the nucleonic origin of the observed TeV emis- 
sion is treated by the community as an inadequate alter- 
native to the IC mechanism, the main argument being the 
low ambient density of the gas in SN 1006, n ~ 0.4 cm^'^ 
(Willingale et al. 1996) as well as the supposed large dis- 
tance to the source of about 2 kpc (Winkler and Long 
Hower these arguments are not sufficiently robust 



1997) 



in order to dismiss such an important possibility with far 
going conclusions concerning the origin of the nucleonic 
component of galactic cosmic rays (Aharonian 19991 ). 

Indeed, the very fact of the existence of 10 TeV elec- 
trons as it follows from the X-ray data is an evidence for a 
strong shock in SN 1006, and implies a large compression 
factor, p ~ 4 or even more, up to 10, as it follows from non- 
linear stu dies o f shock acceleratio n in S NRs ( Malko v 1997, 
Berez hko |l996| , Berezhko & Volk |l997l Volk |l997t Baring 



et al. 1999). Therefore it is not unrealistic if one assumes 
a significantly higher density of the gas in the rim region, 
e.g. n ~ 2cm~'^. The current estimates of the distance to 
the source also contain large uncertainties, and do not ex- 
clude a distance of about 1 kpc or even less. In particular, 
one of the recent studies of SN 1006 based on the ROSAT 
observations suggests a distance d = 0.7 ± 0.1 kpc (Will- 
ingale et al. 1996). To explain the observed flux of TeV 
emission, for Fp = 2 the scaling factor A ~ 0.83 is needed 
(Fig. ^. This requires for the total energy in accelerated 
protons 



p _ 4 X 10^^ (n/2 cm-3)-i dl^, erg . 



(20) 



This is only ~ 10 per cent of the total kinetic energy 
of explosion estimated for SN 1006 as 5 x 10^° erg. The 
numerical calculations within the kinetic theory of shock 
acceleration show (Berezhko |199S ) that even for a large 
distance to SN 1006, the energy budget in accelerated pro- 
tons can be sufficient to explain the TeV 7-ray emission 
from SN 1006 by hadronic interactions. 



10 



o 10 



10 



: 1 

r = 


1 1 
2.1 


1 : 

- 


-1- 






1 / 


2 ^ ^ ' 


A\ : 

\\ " 




1 .8 


V' 


/ / 

- ' / 

- '/ ^ 




V 
% 


/'/ 1 / 


1 1 


i 

1 \ 



9 10 11 
Log(E/eV) 



12 



13 14 



Fig. 7. The fluxes of 7r°-decay 7-rays calculated for the 
proton spectrum in Eq.(l) with E — 200 TeV, and 3 dif- 
ferent power-law indices: Fp = 1.8 (dashed), 2 (solid), and 
2.1 (dot-dashed). The fluxes are normalized to the ob- 
served flux at 3 TeV. This results in the values of scaling 
parameter A = 0.44, 0.83, and 1.35, respectively. 



In fact this estimate of the scaling factor A derived 
from the comparison of calculated and observed TeV 
fluxes depends significantly on the spectrum of protons. 
In Fig. 7 we show the 7r°-decay 7-ray fluxes calculated for 
3 spectra of protons taken in the form of Eq. (1) with spec- 
tral index Fp =1.8, 2, and 2.1. The fluxes are normalized 
to the observed flux at 3 TeV. The corresponding values of 
the scaling factor are A=0.44, 0.83, and 1.35, respectively. 
Even in the case of relatively soft spectrum of protons 
with Fp = 2.1, the required scaling factor is still accept- 
able taking into account that more than 20% of the energy 
of the supernova explosion can be transformed into accel- 



erated protons (Berezhko & Volk 1997). Steeper spectra 
of protons with a slope Fp > 2.1 do not match the energy 
budget of the source. Such spectra are excluded also by 
the EGRET flux upper limit, F-y < 1.7 • lO-'^ph/cm^s at 



E > 100 MeV (Mori |1997| ). 

The 7-ray spectra shown in Fig. 7 are calculated as- 
suming for the maximum energy of protons Eq = 200 TeV. 
Although the exact value of Eq do not noticeably change 
the requirements to the scaling factor A, it has a notice- 
able impact on the spectral form of 7-rays above 10 TeV. 
Therefore only future precise spectroscopic measurements 
in this energy region could provide an important informa- 
tion about Eq. At the same time the existing sub- 10 TeV 
data already tell us that within the framework of nucle- 
onic model of TeV radiation of SN 1006 the cutoff energy 
Eq could not be significantly less than 100 TeV. 

Although in this paper we do not attempt to discuss 
the aspects related to the theory of diffusive shock accel- 
eration in SNRs, it is clear that such large values of Eq 
could be achieved only under the assumption of a large 
magnetic field in the shock region, B ~ 10"'' G or even 
more (see Eq. 11). Obviously this assumption does not 
leave any chance for the IC mechanism to give a notice- 
able contribution into the observed TeV emission. On the 
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Fig. 8. The nonthermal radiation of SN 1006 calculated 
for large magnetic field in the NE rim B = 70 /iG. Solid 
curves correspond to the synchrotron radiation of directly 
accelerated electrons, while the dotted curve corresponds 
to the synchrotron radiation produced by secondary (tt^- 
decay) electrons. The dashed curve shows the flux of tt"- 
decay 7-rays calculated for proton spectrum with Fp = 2, 
Eo = 200 TeV, and scahng parameter A = 0.83. The dot- 
dashed line shows the flux of IC 7-rays contributed mostly 
by directly accelerated electrons. All fluxes shown are well 
dominated by the emission produced in the rim. 

other hand the featureless X-ray radiation could be natu- 
rally explained by the synchrotron radiation of accelerated 
electrons. In such strong magnetic fields the maximum en- 
ergy of accelerated electrons is determined by the balance 
between the acceleration rate and the synchrotron energy 
loss rate: 



Eo ~ 45 7?" 



-1/2 



B 



10 mG 



-1/2 



2000 km/s 



TeV , (21) 



and therefore, as it follows from Eq.(8) the energy of ex- 
ponential cutoff in the synchrotron spectrum. 



Er, 



0.5(i;,/2000 km/sfr]-^ keV . 



(22) 



does not depend on neither the magnetic field nor the age 
of a SNR. The very fact of detection of nonthermal X-rays 
above 1 keV testifies that the particle acceleration takes 
place in the limit of Bohm diffusion (jj 1). 

Therefore an observation of flat X-ray spectrum in the 
region 1-10 keV from any SNR would imply, for a reason- 
able shock speed in the Sedov phase Vg < 4000 km/s, a 
more effective acceleration than it follows from the pre- 
diction of conventional parallel shock acceleration model. 

In the regime of non-negligible synchrotron losses we 
should expect, in addition to the exponential cutoff at 
-Em, another spectral feature ~ a spectral break (increase 
of the power-law index by 0.5) at the energy E^ which is 



determined from the condition t 



to: 



In Fig. ^ we show the spectrum of synchrotron radia- 
tion calculated for the following parameters: B = 70 fiG, 
to — 700 yr, Vg = 2900 km/s, and rj — 1. For these pa- 
rameters the spectral break is at iS'b — 0.02 keV, while the 
exponential cutoff £"111 — 1 keV. Although there is a gap 
by a factor of 20 between these two characteristic pho- 
ton energies, the additional gradual steepening due to the 
exponential cutoff in the electron spectrum becomes no- 
ticeable already at energies E^ > 0.1i?m- 

5. Conclusions 

In the previous theoretical studies of the IC 7-radiation 
of SN 1006 a homogeneous one -zone mode l has been im- 
plied (Mastichiadis 



1996 



Yoshida 



de Jager |l996t Pohl 
& Yanagita 1997 ) which assumes that all accelerated elec- 
trons are conflned in the NE rim. The magnetic field in 
the rim derived from these studies which would not con- 
tradict the observed X-ray to TeV 7-ray fiux ratio should 
not exceed 7 ^G. 

Our study of the 7-ray production in SN 1006 which 
includes the effect of convective and diffusive escape of 
electrons, allows some increase of the magnetic field up 
to 10 fiG, with preferable values being in a rather nar- 
row range between 6 and 8 /iG. Given the short time of 
particle acceleration available, to < isNiooe — 10^ yr, the 
assumption of such a small magnetic field requires a very 
high shock speed in SN 1006, Vg > 3000 km/s, in order to 
provide acceleration of electrons to energies Eo ~ 30 TeV 
needed for explanation of X-ray and TeV observations. 
This implies a large distance to the source, d > 1.5 kpc. 

This model allows definite predictions which could be 
tested by future observations. A significant, if not domi- 
nant, fraction of the IC TeV emission is to be produced 
outside of the NE rim. Namely, a 7-ray flux compara- 
ble with the flux from NE rim, should be then expected 
from the extended inner region of the remnant adjacent 
to the rim. The size of this region depends on the charac- 
ter of propagation of electrons in the rim. Indeed, even in 
the case of diffusion of electrons in the Bohm limit in the 
remnant with magnetic field B2 ^ 2 /iG, the characteris- 
tic distance of penetration of electrons towards the central 
region of the remnant could be estimated as 

0.75 (i^o/lOTeV)^/^^ 



liE,) ~ y/2D{E,)to 



pc. 



(24) 



Eu 



2.9(5/10 ^iG)-'' (to/lO-" yr)"^ keV. 



(23) 



Since the IC production of TeV 7-rays (on 2.7 K MBR) 
with energies less than several TeV takes place in the 
Thompson regime, and therefore their characteristic en- 
ergy scales as S ~ 1 (i;e/20TeV)2TeV, the size of the 
emission region of TeV gamma rays should be as large as 
the width of the rim, and increases linearly with energy 
of 7-rays. Moreover, if the propagation of electrons in the 
remnant is much faster than in the Bohm limit (ji ^ 1), 
the electrons could fill up practically the entire remnant. 
In that case the energy dependence of the size of the 7-ray 
emission will be significantly weakened. 
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Another distinctive feature of the IC origin of TeV ra- 
diation is the spectral shape of the radiation: very hard, 
with the photon index ~ 1.5 below 1 TeV, a flat spec- 
trum with ~ 2 between 1 and 10 TeV, and very steep 
above 10 TeV. 

The alternative to the IC interpretation is the nucle- 
onic origin of the observed TeV radiation. This interpreta- 
tion becomes energetically comfortable, which implies no 
more than 10^° erg in accelerated protons and nuclei, if we 
assume (i) a small distance to the source of about 1 kpc or 
less, and (ii) a significant enhancement of the gas density 
in the rim due to strong shock compression. This predicts 
a compact size of the TeV emission comparable with that 
of the X-ray rim. In the region below 1 TeV the spec- 
trum of 7-rays of nucleonic origin, with a power-law index 

~ r ~ 2 is expected steeper than the spectrum of IC 
7-rays. However a flatter spectrum of 7r°-decay 7-rays can- 
not be excluded if the protons would ha ve an acceleration 
spectrum as hard as E~^-^ (see Malkov 1997 ). 

The shape of the spectrum of 7r°-decay 7-rays above 
1 TeV depends on the characteristic maximum energy of 
accelerated protons. Since the flux of TeV 7-rays is no 
more connected with the X-ray fluxes, the value of the 
magnetic field in the acceleration region (NE rim) could 
be adopted as high as ~ 100 /xG. Correspondingly, the 
energies Eq > 100 TeV could be achieved. On the other 
hand if Eq would be significantly less than 100 TeV, one 
may expect a turnover in the spectrum of 7-rays above 
several TeV. 

Thus, it is quite possible that both in the low (sub- 
TeV) and in the high (multi-TeV) energy regions the IC 
and TT^'-decay 7-rays may have similar spectra. Therefore 
the spatial rather than spectroscopic measurements of 7- 
radiation above 100 GeV with future lACT arrays could 
provide decisive information about the origin of very high 
energy radiation of SN 1006. The IC models predict sig- 
nificant 7-ray fluxes not only from the rim, but also from 
the inner parts of the remnant. Meanwhile, the 7r°-decay 
7-rays trace the density profile of the gas in the produc- 
tion region, therefore one should expect a rather compact 
7-ray source essentially coinciding with the rim. 

We emphasize that both the electronic and the nucle- 
onic models of the observed TeV radiation of SN 1006 re- 
quire a very effective acceleration of particles in the regime 
close to the Bohm limit. Also, the both models predict a 
very high conversion efficiency of the total kinetic energy 
of the explosion available (~ 5 x 10^*^ erg) into relativis- 
tic particles. The energy required for relativistic electrons 
within the framework of the IC model is predicted with 
a good accuracy, given the small range of allowed B-field: 
Wc ~ 5 X 10"*^ erg. This is by two orders of magnitude 
larger than the energy contained in the magnetic field; 
this suggests that the conditions in SN 1006 are far from 
the equipartition regime. 

The nucleonic model of TeV emission requires total 
energy in accelerated protons from 2 x 10^" erg down to 



10 erg, depending on the enhancement factor for the gas 
density in the rim, the spectral index of the accelerated 
protons, and the distance to the source. In particular for 
moderate assumption for the spectral index Fp = 2, and 
the gas compression ratio p = 4 (i.e. n = 1.6 cm~'^), the 
energy in protons is estimated « 5 x 10'^®(d/l kpc)^ erg. 
Note that since the magnetic field in the rim for the nucle- 
onic model could be as high as 100 fxG, this model predicts 
a situation close to the equipartition. 

High magnetic fields in the rim can equally well explain 
the current X-ray data. The synchrotron X-radiation in 
this case is produced in the regime of saturated energy 
losses of electrons, which results in a somewhat different 
spectral shape of the synchrotron fluxes. In both large 
and small magnetic field limits, the spectrum would be 
best fitted in the form F{E) oc i;-"- exp[-(£;/£'*)-i/2], 
where ~ 2.1 and ax — 1-6 for the cases of large and 
small magnetic fields, respectively. 
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